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Abstract 

We describe our theoretical formalism and computational scheme for mak- 
ing ab-initio calculations of the dynamic paramagnetic spin susceptibilities of 
metals and alloys at finite temperatures. Its basis is Time-Dependent Den- 
sity Functional Theory within an electronic multiple scattering, imaginary 
time Green function formalism. Results receive a natural interpretation in 
terms of overdamped oscillator systems making them suitable for incorpora- 
tion into spin fluctuation theories. For illustration we apply our method to the 
nearly ferromagnetic metal Pd and the nearly antiferromagnetic chromium al- 
loy CV95V5. We compare and contrast the spin dynamics of these two metals 
and in each case identify those fluctuations with relaxation times much longer 

than typical electronic 'hopping times'. 
75.40Gb,75.10Lp,71.15Mb,75.20En 
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There is renewed interest in spin fluctuations in materials close to magnetic order. This 
is due in part to a realisation that nearly critical magnetic fluctuations may be important 
factors governing the non-conventional properties of a wide range of materials which in- 
clude the high T c superconducting cuprates and heavy fermion systems JIJ. The strongly 
correlated electrons in many of these systems, however, have meant that most theoretical 
work has concentrated on parameterised models in which the electronic motion is treated 
rather simply. Another complementary approach is to use an ab initio theory such as Time 
Dependent Density Functional Theory (TDDFT) 0] but apply it to materials where it can 
be expected to work i.e. where the effects of electron correlations are not so important, 
but which otherwise have important similarities to the systems in question. For example, 
with its perovskite structure containing transition metal(TM)-oxygen planes, Sr 2 RuOi has 
several aspects in common with the HTC materials. But the presence of the 4d TM Ru 
rather than the narrower band 3d TM Cu means that electron correlation effects are smaller 
and therefore DFT-based calculations can provide a valuable starting point. Moreover, its 
exotic superconductivity at low temperature seems likely to be affected by spin fluctuations 
0-0. Concerning another example, the transition temperature separating paramagnetic 
and magnetically ordered phases of the cubic transition metal compound MnSi, which has 
the B20 crystal structure, is driven down to zero temperature upon the application of pres- 
sure. In the vicinity of the critical pressure for this quantum phase transition the system 
exhibits non- or marginal Fermi liquid properties . In this paper we investigate the energy 
and temperature dependence of spin fluctuations in two systems which although structurally 
and compositionally simpler than the two mentioned above, are also close to magnetic phase 
transitions. The first is palladium which is nearly ferromagnetic and the second is the nearly 
anti-ferromagnetic CV95V5 alloy. 

Theoretical models in which an effective action for the slow spin fluctuations is written 
down have contributed greatly to our understanding of the properties of itinerant electron 
systems close to magnetic order 0,0]. Recent work which can incorporate results from 
DFT-based 'Fixed Spin Moment' (FSM) electronic structure calculations, treats these fluc- 
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tuations classically ||. A Landau- Ginzburg-like energy functional is written down, a free 
energy constructed which includes terms describing the interactions between the fluctuations 
(the mode-mode coupling) and properties such as the static susceptibility, specific heat and 
resistivity calculated. The FSM electronic structure calculations can be used to determine 
the coefficients in this functional ||. Many informative studies have been carried out. For 
those calculations with such a DFT basis, these still remain qualitative investigations be- 
cause of the lack of a prescription for the effective number of modes to include in the theory 
and its variation with temperature Stoner single-particle excitation effects are also 

largely ignored. Both these issues can be addressed by the development and application of 
methods to calculate the temperature-dependent, dynamic paramagnetic spin susceptibility 
of nearly magnetic materials. The development of dynamic susceptibility calculations is 
particularly pertinent now that inelastic neutron scattering experiments, such as the time of 
flight measurements, have developed to the extent that spin fluctuations in nearly magnetic 
metals can be accurately measured ||. 

We have recently devised and proven a new scheme for calculating the wave-vector and 



frequency dependent dynamic spin susceptibility of metallic systems [ID which is based on 
the Time Dependent Density Functional Theory (TDDFT) of Gross et al. and as such is 
an all electron theory. This enables the temperature dependent dynamic spin susceptibility 
of metals and compositionally disordered alloys to be calculated. 

Over the past few years great progress has been made in establishing TDDFT M. Analogs 



of the Hohenberg-Kohn [11 theorems of the static density functional (DFT) formalism have 
been proved and rigorous properties found. By considering a paramagnetic metal subjected 
to a small, time-dependent external magnetic field, b(r, t) which induces a magnetisation 



m(r, t) and using TDDFT in |10[ an expression for the dynamic paramagnetic spin sus- 



ceptibility x(.Q.i w ) via a variational linear response approach can be derived [O]. Accurate 



calculations of dynamic susceptibilities from this basis have been scarce (e.g. (13]]) because 
they are difficult and computationally demanding. In ref. |]TU[ we showed that these prob- 
lems can be mitigated by accessing x(q, w) via the corresponding temperature susceptibility 



x(q, w n ) where w n denotes a bosonic Matsubara frequency |p^4]| . We demonstrated this ap 



proach by an investigation into the nature of the spin fluctuations in paramagnetic Cr and 
compositionally disordered CV95V5 and Crg 5 Re^ alloys with good agreement with experi- 
mental data |P-0|| . For example, recent inelastic neutron scattering experiments ||[15| have 



measured incommensurate AF 'paramagnons', persisting up to high frequencies in the latter 
system which were also shown in our calculations 



In [|10| we provided a brief outline of this approach for crystalline systems with one atom 
per unit cell only so in the next section we provide its generalisation to multi-atom per unit 
cell materials in which one or more sublattices may be substitutionally disordered. We also 
describe the techniques involved in the method in more detail. In the following section we 
discuss our study of the frequency, wave-vector and temperature dependence of the spin 
fluctuations in Pd and identify the paramagnon regime. We then move on to carry out a 
similar study of the incommensurate antiferromagnetic spin fluctuations in the dilute CV95V5 
alloy. In each case the results can be interpreted simply in terms of an over-damped harmonic 
oscillator model 0. In the latter case we show how the parameters scale with temperature 
and compare and contrast the two systems. In the final section we draw some conclusions 
and remark how this work has the potential to be integrated into spin fluctuational theories 
of nearly and weakly magnetic metals ||. 

I. THE DYNAMICAL SPIN SUSCEPTIBILITY x(Q,W). 

The equilibrium state of a paramagnetic metal, described by standard DFT, has density 
Po(r) and its magnetic response function 

x(rt,rt) - WJ) 1 ,po ( } 

is given by the following Dyson-type equation. 
X (rt;r't') = X s(rt;r't') + J dr, J dt x J dr 2 J dt2X.{rtiT 1 t 1 )K m (r 1 t 1 -,r 2 t 2 )x(r 2 t 2 ,r f t') (1.2) 
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Xs is the magnetic response function of the Kohn-Sham non-interacting system with the 
same unperturbed density p as the full interacting electron system, and 

is the functional derivative of the effective exchange- correlation magnetic field with respect 
to the induced magnetisation. As emphasised in ref. || eq.1.2 represents an exact represen- 
tation of the linear magnetic response. The corresponding development for systems at finite 
temperature in thermal equilibrium has also been described ||12|| . In practice approximations 
to K xc must be made and this work employs the adiabatic local approximation (ALDA) |2j 
so that 

K xc rt) = ^ lp 0i m=o<3(r -T)6(t-t) 

= I xc (r)6{r-r')6(t-t') (1.4) 

On taking the Fourier transform of eq.1.2 with respect to time we obtain the dynamic spin 
susceptibility x( r , r '; w). 

For computational expediency we consider the corresponding temperature susceptibil- 
ity |14| x(r, r'; w n ) which occurs in the Fourier representation of the temperature function 
x(rr; rV) that depends on imaginary time variables t,t' and w n are the bosonic Matsubara 
frequencies w n = 2n7rfcgT. Now x(r,r';w n ) = %(r,r';zu; n ) and an analytical continuation 
to the upper side of the real w axis produces the dynamic susceptibility x( r ; r '; w). 

We define our general system as having a crystal structure with lattice vectors {Ri} and 
where there are N s non-equivalent atoms per unit cell. On the I th of the N s sublattices there 
are Ni atomic species with concentrations (an = 1 , - - - , iVj) . In each unit cell the N s 
atoms are situated at locations aj, Z = 1, • • • , N s . The volume of the unit cell is Vws- On 
carrying out a lattice Fourier transform over the lattice vectors {Rj} we obtain the following 
Dyson equation for the temperature susceptibility 

Nl ' 

X ai (x,,x;,,q,u; n ) = J2 ^XT 11 ' xj,, q, w n ) 

7;' 
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with x;,Xj, and xj« measured relative to atomic cells centred on a z , ay and a;// respectively. 
Xs llv describes the non-interacting susceptibility averaged over all configurations in which 
one site on sublattice I is occupied by an atom and another site on sublattice V is occupied 
by a 7;/ atom. is the full susceptibility averaged over all configurations where an ai atom 
is located on one site on sublattice I. 

In terms of the lattice Fourier transform of the DFT Kohn-Sham Green function of the 
static unperturbed system the non-interacting susceptibility can be written approximately 

m 



X? 7i '(x«,x;,,q,u; n ) = 
1 f dk 

-75 Y. / (C(x z , x,,, k, fi + iu m )) aa (G(x[,, xj, k + q, /i + + w n ))) % , ai (1.6) 

P m V BZ 

where the integral is over the Brillouin zone of the lattice and k, q and k + q are all 
wavevectors within this Brillouin zone which has volume vbz- is the chemical potential 
and fermionic Matsubara frequency (2n + l)nkBT. 

The Green function can be obtained within the framework of multiple scattering (KKR) 
theory |L6| and this makes this formalism applicable to disordered alloys as well as ordered 
compounds and elemental metals, the disorder being treated by the Coherent Potential 
Approximation (CPA) |T7|]. Then the partially averaged Green function with species ol\ on a 
site at Rj + sq on the I th sublattice and species on a site at Rj + on sublattice number 
/' 

(G(Ri + 3li+ x,, Rj + 3l v + x r , z)) aitll , (1.7) 
can be evaluated in terms of deviations from the Green function of an electron propagating 



through a lattice of identical potentials determined by the CPA ansatz [18]. The lattice 
Fourier transform used in equation ( |1.6| ) is expressed as 

(G(Xj,3^,,k,z)) aj7j< = J2 Z L^hZ)(Tl l ' L ,(\L,z)) ani ,Z^{x! v ,z) 

L,L' 
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- Zt\^z){r^ LI {z)) ani Zl^ u z^ (1.8) 

where Z^ 1 and J^ L are respectively the regular and irregular solutions of the Schrodinger 
equation in an atomic cell on sublattice I inhabited by an on atom and L,L' represent 
the angular momentum quantum numbers. The lattice Fourier transform of the averaged 



scattering path operator r l £ l L ,(Ri — Rj, z) is specified by the following expressions [RJ 



{T l £ l tL ,(k,z)) ani , = £ D l auLLi (z)f^ L2 (k,z)D i ; iiML , (1.9) 
L1L2 

where, in matrix notation with respect to angular momentum L,(D describes the transpose) 

D l auLLl (z) = [l+f u (z)(m^(z)-m c ' l (z))]ll l (1.10) 

with m ai being the inverse of the scattering t-matrix t ai of an oli atom potential, m c ' 1 the 
inverse of the CPA t-matrix for the I th sublattice, t c ' 1 and f 11 ' the unit cell-diagonal part of 
the CPA scattering path operator. The following averages also involve these quantities 

(^'v{z)) aail = £ ^M^mMKmuW 

L1L2 
Li 

Finally the lattice Fourier transform of this CPA scattering path operator is given by 

< L ,(k, z) = tH L ,(z)8 w + E E ti Ll {z)gZ L ^ z)fZ>(K z) (1.12) 

i" L1L2 



and g l £ lL2 (k, z) are the structure constants for the crystal structure |19| . 



To solve equation ( |1.5| ), we use a direct method of matrix inversion in which, for example, 
Xs is cast into matrix form of order (J2iti S1N1) x (J2iti S1N1) where Si is the number of spatial 
grid points associated with the I th sublattice. Local field effects are thus fully incorporated. 
The full Fourier transform 



X 



[q,q;w n ) = (l/VW)EE<E eMHa! ~ ai,) / d * I d^ v ^-*'i<h ai M,,^Wn) (1.13) 

l a t V 



can then be constructed. 

The most computationally demanding parts of the calculation are the convolution inte- 
grals over the Brillouin Zone which result from the expression for Xs, eq. (|1.6|) i.e. 



Jdkf l l' Li (k,z)fl l 2L3 (k + c i ,z') (1.14) 

Since all electronic structure quantities are evaluated at complex energies z, these con- 
volution integrals have no sharp structure and can be evaluated straightforwardly by an 
application of the adaptive grid method of B.Ginatempo and E.Bruno [^] which has been 
found to be very efficient and accurate. In this method one can preset the level of accuracy 
of the integration by supplying an error parameter e. In the calculations described in this 
paper we have used e = 10 -4 so that we obtain 4 significant figure accuracy on Xs- This 
control on accuracy is crucial for the description of the long wavelength paramagnons in Pd. 
In this case we needed to sample around 2000 k-points at energies on the contour distant 
0.5Ry. off the real axis but up to 50000 k-points very close (0.001 Ry.) to the real energy 
axis. Around half that number of points were required for the CV95V5 calculations. 

We evaluate the Matsubara sum shown in equation ( |1.6| ) by using suitable complex 
energy contours to enclose some of the Matsubara poles and exploiting the hermiticity of 
the single-electron Green function G(r, r', z). Consequently a quantity like Xs at a positive 
bosonic Matsubara frequency w n with form 

A(r,r',w n ) = -^-^2G(r,r',fi + iu m )G(r',r,fi + iu m + iw n ) (1.15) 

r m 

where the fermionic Matsubara sum is over m = 0, ±1, ±2, ■ • • , ±00 can be rewritten as 

2 Mm ax 

A(r, r', w n ) = — - Re(G(r,r',fi + iu m )G(r',r,fi + iu m + iw n )) + (1.16) 

P m=0 

If 1 n-1 

-Im / dz(G(r, r',z)G(r',r, z + iw n ))f(z) - -V G(r, r', \i - iv m )G(r', r, \i + iw n - iv m ) 

Both sums are over positive fermionic frequencies only and f(z) is the Fermi function. C is 
a contour which encloses the first M max fermionic Matsubara poles with the real axis. In our 
calculations we have used a simple box contour whose top lies exactly halfway between two 
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neighbouring i/ m 's which makes the Fermi factor at this z real. The terms in the Matsubara 
sums for z/ m 's > 0.001 are evaluated by interpolating from a grid v = 0.001 x 10 aip , with 
p — 0, 1, • • • , 30 or are evaluated directly. M max is chosen so VM max ~ 1. 

Once the temperature susceptibility x(q, q; w n ) has been calculated the dynamic sus- 



ceptibility can be found. As discussed in ref. [|14||, for example, we can define the retarded 



response function x(q, q, z) of a complex variable z. Since it can be shown |14| formally 
that lini^oo x{ z ) ~ 1/z 2 and we can obtain x(iu> n ) from the above analysis it is possible 
to continue analytically to values of z just above the real axis, i.e. z = w + irj. In order to 
achieve this we fit our data to a function describing a set of overdamped oscillators, i.e., 

X(q, q, w n ) = 2^ j— —, 1.17 

y (1 + (w n /r M (q)) + (w n /fi M (q)) 2 ) 

in which M is an integer. From this the dynamic susceptibility can be written down and 
the imaginary part becomes 

m It 1 - (^AMq)rr + (^/ r M(q)) 2 } 

This form ensures that the sum rule involving the static susceptibility x(q) is satisfied, i.e. 

2 f°° Jmx(q,q,w) 

X(q) = - / dw 1.19 

n Jo w 

In nearly magnetic metals the nature of the spin fluctuations can be succinctly described 
via the variance < m 2 (q) >. From the fluctuation dissipation theorem 

/ 2/ w 1 f°° A Jm x(q, q,w) onA 

< m (q) >= - / dw—- ^r— (1.20) 

v ; W-oo (l-e-^) v y 

By introducing an energy cutoff in this expression the slower fluctuations can be revealed. 
We find that very good fits are obtained with small M in equation (|1.17| ) for a wide 



range of w n J s ( < 500 meV) for the systems we have examined to date. For Cr and its dilute 
alloys using M — 1 in equation ( |1.17| ) provides an excellent fit. This means that a close 



analogy with an overdamped harmonic oscillator can be made. In the case of palladium, 
equation (|1.17| ) with M = 2 gives an excellent representation of x(q, q, u; n ). Now the 



spin fluctuations can be described in terms of two overdamped oscillators with two distinct 
timescales. The slower one describes the famous paramagnons of palladium. 
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II. THE SPIN FLUCTUATIONS OF NEARLY FERROMAGNETIC PALLADIUM. 



Although possessing one of the highest Fermi energy densities of states of all the tran- 
sition metals, palladium has neither a magnetically ordered nor a superconducting ground 
state. With its nearly filled d-bands, it is, however, on the verge of being a ferromagnet. 
This is shown by its uniform paramagnetic spin susceptibility which is greatly enhanced 
over that of an equivalent non-interacting system by exchange-correlation effects. Moreover, 
the temperature dependence of the susceptibility also demonstrates effects of long-lived fer- 
romagnetic spin-fluctuations or paramagnons. It has even been proposed that these spin 
fluctuations might induce p-wave superconductivity in this material on account of similari- 
ties with spin-fluctuation induced superfluidity in He 3 [21 . 



Several theoretical studies of the spin fluctuations of Pd have been carried out. The early 
ones were based on a jellium model for the electrons and examined the contribution that 
spin fluctuations make to the specific heat and their influence on any superconducting phase 
transition as well as describing the cross-section that would be measured in an inelastic 
neutron scattering experiment |22]. Later work has attempted to remove the deficiencies 
of such a single-band model and to describe the band structure effects more realistically 
||13||. From his calculations of the inelastic neutron scattering cross-section for a range of 



scattering angles and incident energies, Doniach |S2[ estimated the circumstances under 
which the paramagnons of Pd might be observed. Very recently time of flight measurements 
on a neutron spallation source have succeeded in detecting these long-wavelength, long-lived 
spin fluctuations |§. 

Using the formalism described in the last section, we have calculated the temperature 
susceptibility of Pd at T = 100K and have then determined the dynamic susceptibility 
x(q, w). Currently our calculations use descriptions of the equilibrium states in terms of 
scalar-relativistic atomic sphere approximation (ASA) effective one-electron scattering po- 
tentials and charge densities. The formalism, however, is applicable to full potentials and 
non-spherically symmetric charge densities. The lattice spacing of 3.83 A which we used in 
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the Pd calculations was determined from the minimum of the ASA total energy and found 
to be 1.5% smaller than the experimentally measured value. At 15.7 fi^eV^ 1 our calculated 
uniform susceptibility x(0, 0) is somewhat lower than the value 24.3 [i 2 B eV~ l measured ex- 
perimentally ||. For small wavevectors q the inverse of our calculated static susceptibility 
depends on the magnitude squared of q, i.e. 0) ~ X _1 (0, 0) + cq 2 as also found exper- 

imentally and in other calculations ||13|| . We obtain a value of approximately 900 [i 2 B AmeV 



for c whereas the value extracted from inelastic neutron scattering experiments is 294 ± 130 
fi 2 B AmeV. Since palladium is so close to a ferromagnetic phase transition both x(0, 0) and 
c are very sensitive to the value of the lattice parameter used and the RPA approximation 
used for the exchange-correlation interactions. Using the experimental lattice spacing in the 
calculations greatly increases x(0, 0) (to 23.5 jji 2 B eV~ l close to the value determined from 
experiment) and both x(0, 0) and c are likely to be altered if spin fluctuation interaction 
effects are treated beyond the LDA [|7| . Key aspects of the dynamic behaviour turn out not 
to be so sensitive. 

In figure 1 we show the calculated imaginary part of x(q, w) for a range of frequencies up 
to 200 meV for selected wavevectors q along the (1, 0, 0) direction. The sharp peaks at low 
frequencies for q < 0.1 indicate paramagnon-like behaviour. This is illustrated further in 
figure 2 which shows the magnetic correlations < m 2 (q) > using energy cutoff's of 500meV 
and 50meV. Once again the paramagnons are clearly visible for a narrow region of small 
wave-vectors q. 

We find that the results can be interpreted in terms of two overdamped simple harmonic 
oscillators each with a characteristic timescale. The longer timescale one encapsulates the 
paramagnon behaviour. Jmx(q, q, w) can now be written in the following way 

(1 + (w/r^q)) 2 ) (1 + (w/r 2 (q)) 2 ) 

Figure 3(a) shows the quantities Xi( t l); X2(q) and their sum whereas figure 3(b) shows 
rf x (q) and T^q). Evidently at small q the slow relaxation process is dominant and 
Ti/Yi(q) describes the long time taken for a spin density perturbation to relax back to the 
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equilibrium value. Moreover, we find that the product of Xii^l) with Ti(q) varies as 7og for 
small q, so that spin density disturbance decays over a time proportional to the wavelength, 
i.e. exhibits Landau damping [|]. We find 70 to be 2.13 fi 2 B A which is in good agreement 
with the experimentally determined value of 1.74 ± 0.8 || and evidently fairly insensitive 
to the treatment of exchange-correlation effects. 



III. INCOMMENSURATE ANTIFERROMAGNETIC SPIN FLUCTUATIONS IN 

C-R95V5. 

Chromium, an early 3d transition metal, loses its anti-ferromagnetic ground state when 
electrons are removed as a suitable dopant is added. For example a strongly exchange- 
enhanced paramagnet exhibiting anti-ferromagnetic paramagnons is formed when just a few 
atomic percent of chromium are substituted by vanadium. In this section we explore the 
similarities and differences between these excitations in one such system, CV95V5, and those 
reported for palladium. We also examine the temperature dependence. 

There is an extensive literature on chromium as the archetypal itinerant anti-ferromagnet 
(AF) and its alloys [fL5|]. It is well known that Cr's famous incommensurate spin density 
wave (SDW) ground state is linked to the nesting wave-vectors q nes t identified in the Fermi 
surface. This feature also carries over to its dilute alloys which have a range of AF proper- 
ties [[31| (see figure 4). Indeed the paramagnetic states of some of these alloys have recently 
acquired a new relevance on account of analogies drawn with the high temperature su- 



perconducting cuprates especially (La c Sri^ c ) 2 Cu04 : |P3| . Starting with 'parent' materials 
CrgsMris or Cr 95 Re^, which are simple commensurate AF materials and then lowering the 
electron concentration by suitable doping causes the Neel temperature, Tjy, to drop rapidly 
to zero. The paramagnetic metal which forms for dopant concentrations slightly in excess 
of the critical concentration for this quantum phase transition is characterised by incom- 



mensurate paramagnetic spin fluctuations. In [10] we described the spin fluctuations in the 



paramagnetic state of three representative systems Cr 95 Re 5 , Cr and Cr 95 V 5 . We tracked 
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the tendency for the spin fluctuations to become more nearly commensurate with increas- 
ing frequency and also with increasing temperature ]T(| as observed in experimental data. 
Our estimates of T/v in Cr^Re^ and Cr of 410K and 280K, respectively, were also in fair 
agreement with the experimental values of 570K and 31 IK. 

Although there have been several simple parameterised models to describe the magnetic 
properties of Cr and its alloys f24| , these have all concentrated on the approximately nested 



electron 'jack' and slightly larger octahedral hole pieces of the Fermi surface (FS) jlEj . These 
sheets can be seen clearly in figure 4 which shows the Fermi surface of CV95V5. The FS of 



Crg^Re^s is close to being perfectly nested and we calculated in [T(| its spin fluctuations 
above T/v to be nearly commensurate for small frequencies w. We also looked at Cr and 
Cr 95 V 5 . Cr's and Cr 95 V^s Fermi surfaces are progressively worse nested than Cr^Re^s. 
In our calculations for these two systems |10[ we found the dominant slow spin fluctuations 
to be incommensurate with wave-vectors equal to the FS nesting vectors q nes t- Cr^V^s is 
shown on figure 4. At best, simple parameterised models only include the effects of all the 
remaining electrons via an electron reservoir. Whilst finding the obvious similarities from 
the FS basis between our results and results from such models, we showed that a complete 
picture is obtained only when an electronic band-filling effect which favours a simple AF 
ordering at low temperature is also considered. 

As for the case of Pd we found that the spin fluctuations in the paramagnetic state of the 
three systems are given an accurate description in terms of a overdamped diffusive simple 
harmonic oscillator model. Here, however, a single channel is sufficient i.e. the susceptibility 
closely fits the following 

Imx(q, q, w) « /f^/^ffi^L 
(1 + (w/r(q)) 2 ) 

In this section we concentrate on the exchange-enhanced paramagnet CV95V5 and focus 
in particular on the temperature dependence of the incommensurate spin fluctuations and 
show how the parameters of the oscillator model scale with temperature. 

Figure 1(a) of reference |Kj showed Jmx(q, q, w) for wavevectors q along the {0,0,1} 
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direction for frequencies up to 500meV. It showed, in agreement with experiment ||, in- 
commensurate AF paramagnons persisting up to high frequencies with intensity comparable 
to that at low w. This is in striking contrast to the paramagnons of Pd. Figure 5 shows 
this persistence of the spin fluctuations in CV95V5 to high frequencies in a concise way. It 
shows the magnetic correlations < m 2 (q) > at 100K for the same two energy cutoffs 50 
and 500meV as was used to obtain figure 2 for palladium. Unlike those in the nearly fer- 
romagnetic metal, the magnetic correlations grow by more than an order of magnitude for 
all the wave-vectors shown as the energy cutoff is extended by a similar factor. The varia- 
tion with temperature of the dynamic susceptibility and the magnetic correlations have also 
been investigated. Figure 6 shows < m 2 (q) > calculated with a energy cutoff of 500meV 
for 50K, 300K and 600K. There is a general trend for the weight to transfer towards the 
commensurate wavevector (0, 0, 1) with increasing temperature. For a smaller energy cutoff 
< m 2 (q) >, integrated over all q, decreases with increasing temperature. 

Figures 7(a) and (b) show that both the static susceptibility x(q) and the spin relaxation 
time r _1 (q)/ft also peak at the incommensurate nesting vector (0,0,0.9) at low tempera- 
ture and that weight is shifted towards (0, 0, 1) at higher temperatures. Figure 7(b) when 
compared with the analogous figure 3(b) for Pd reveals the main difference between the 
dominant spin fluctuations in the nearly antiferromagnetic CV95V5 and the nearly ferromag- 
netic Pd. For the dominant spin fluctuations, the spin relaxation times are some 50 times 
faster. The small Q behavior is different also. (Q = q — qo, where qo = for Pd and an 
incommensurate wave-vector for CV95V5.) The product 7(q) of x{°d an d T(q) tends to a 
constant for values of q close to qo for CV95V5 (consistent with ideas that the dynamical 
critical exponent should be 2 for antiferromagnetic itinerant electron systems) unlike the 
linear variation of 7(q) with |q| for small q found in Pd. Between 50K and 300K, 7(q) 
shows a weak variation with temperature where it has values of roughly 4 fi B and 6.5 fi 2 B at 
the incommensurate wave-vector q and q = (0, 0, 1) respectively. Above 300K, 7(q) loses 
the trough around qo and by 600K it is roughly constant at 5.5 fi B for the range of q from 
qo to (0,0,1). 
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IV. CONCLUSIONS 



We have described the framework and method of implementation of our scheme for car- 
rying out ab-initio calculations of the dynamic paramagnetic spin susceptibilities of solids at 
finite temperatures. The approach is based upon Time-Dependent Spin Density Functional 
Theory and is applicable to compositionally disordered alloys with multi-atom per unit cell 
crystal structures. (We note here that this approach may also be adapted to the study of 
magnetic excitations in magnetically ordered materials.) From an imaginary time multiple 
scattering Green function formalism an expression for the temperature susceptibility has 
been derived and the techniques appropriate for its evaluation described. We have shown 
how the dynamic susceptibility is obtained from this by analytic continuation from Matsub- 
ara frequencies in the complex plane to the real frequency axis ||25|| . This step provides a 
natural interpretation of the spin dynamics in terms of overdamped oscillator models. 

Although ultimately aiming to investigate the spin fluctuations in systems with complex 
lattice structures such as SV^-RwC^ and MnSi which are of topical interest, here we have 
compared and contrasted the nearly ferromagnetic transition metal Pd with the nearly anti- 
ferromagnetic dilute Cr alloy, CV95V5. In both cases we have been able to identify 'slow' 
spin fluctuations so that in due course the mode-mode coupling amongst them may be 
incorporated into spin-fTuctuational theories which describe the low temperature properties 
of these materials HH. For the case of Pd the spin susceptibility can be interpreted in 
terms of two oscillator models. One describes spin fluctuations with fairly 'fast' relaxations 
times but the other describes the spin fluctuations for all wave-vectors which have spin 
relaxation times more than 40 times as slow as the time taken for an itinerant d-electron to 
hop between lattice sites. Clearly there is a natural time scale separation between these spin 
fluctuational modes and the electronic degrees of freedom. Consequently spin fluctuational 
theory can be used with these modes without the need for any wave-vector cutoff. Work 
is in progress to carry this out. On the other hand Cr^V^s dynamic spin susceptibility is 
interpreted in terms of a single oscillator. Here the spin relaxation times are only some 10 
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times slower than typical d-electron hopping times for modes with q's in a limited region of 
the Brillouin zone and therefore mode- mode coupling effects may not be so important. This 
may explain why our calculations |10[ of the static susceptibility of this Cr alloy (and Cr 
and its other dilute alloys) receive an adequate description in terms of what is essentially 
an ab-initio Stoner theory. 
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FIGURES 

FIG. 1. /mx(q, q,w) of Pd at T =100K in units of n B eV~ l for wave-vectors q along the 
{0,0,1} direction where q is in units of n/a (a is the lattice spacing). The frequency axis, w, is 
marked in meV. The full curve is for q =(0,0,0.02), then for curves with decreasing dash-length, 
q= (0,0,0.1), (0,0,0.2) and (0,0,0.5). 

FIG. 2. The variance of the spin fluctuations in Pd at T =100K, < m 2 (q) >, in /i B for 
wave-vectors q along the {0, 0, 1} direction with energy cutoffs of 500 meV (full line) and 50 meV 
(dashed line). 

FIG. 3. (a) The static susceptibility of Pd at T =100K partitioned into 2 parts Xi( ( l)) X2(q) of 
equation ([O]) in units of n B eV~ l for wave-vectors q along the {0,0, 1} direction. xi(q) + X2(q) 
is shown by the full line, Xi{°L) by the long-dashed line and X2(q) by the short-dashed line, (b) 
The time scales from eq. ( |2.1| ) (divided by Ti) in eV -1 , rj~ 1 (q) (full line), r^" 1 (q) (dashed line) for 
Pd at T =100K. Note T^^q) has been scaled up by 100 for display purposes. 

FIG. 4. The Fermi surface (FS) of CV95V5 from Bloch spectral function calculations in the 
Qz = plane. The dashed arrows mark the nesting vector q raes t connecting the -ff-centered hole 
octahedron to the slightly smaller T-centered electron surface. The FS is well defined with only a 
small disorder broadening < 0.0l7r/a. 

FIG. 5. The variance of the spin fluctuations in CV95V5 & t T =100K, < m 2 (q) >, in fi B for 
wave-vectors q along the {0,0, 1} direction with energy cutoffs of 500 meV (full line) and 50 meV 
(dashed line). 

FIG. 6. The variance of the spin fluctuations in CV95V5, < m 2 (q) >, at T =50K (full line), 
300K (long-dashed) and 600K (short-dashed), in [i 2 B for wave-vectors q along the {0, 0, 1} direction 
with a energy cutoff of 500 meV. 
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FIG. 7. (a) The static susceptibility of Cr 95 V 5 at T =50K (full line), 300K (long-dashed) and 
600K (short-dashed), in ^eV -1 along the {0,0, 1} direction, (b) The time scale (divided by fi) in 
eV~ l , T _1 (q) for CV95V5 at the same temperatures. 
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